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We develop an efficient method based on the linear regression algorithm to probe the cosmological 
CPT violation using the CMB polarisation data. We validate this method using simulated CMB 
data and apply it to recent CMB observations. We find that a combined data sample of BICEP 1 
and BOOMERanG 2003 favours a nonzero isotropic rotation angle at 2.3<r confidence level, i.e., 
a = —3.3° ± 1.4° (68% CL) with systematics included. 

PACS numbers: 


I. INTRODUCTION 

The Charge-Parity-Time (CPT) Reversal invariance is a fundamental symmetry in the standard model (SM) of 
particle physics, so a CPT-violating signal is a smoking gun of new physics beyond SM. A number of ground-based 
experiments have been built in order to test this symmetry to high precision. However, no statistically significant 
signals have been found from these lab-scale experiments so far. 

On the other hand, the CPT symmetry needs to be reexamined on cosmic scales. In the universe, the CPT symmetry 
can break down in the model with an effective coupling £j nt = between a non-conserved current J M of the 

SM particles and a scalar field q !>, which may be identified as the dark energy or the Ricci scalar R , and / is an 
arbitrary function of (j). As an external field, the dynamics of <fi is fixed, therefore the Lorentz and CPT symmetries of 
the SM particles break down spontaneously. This has been used to interpret the matter-antimatter asymmetry with 
the current being the baryon or lepton currents [1-3] . In these models the CPT-violating signal is large at the early 
stage to generate enough baryon number asymmetry but too small to be detectable by the ground based experiments 
at present time. 

Another interesting phenomenon of the cosmological CPT violation arises if the is the Chern-Simons current 
of the electromagnetic field. In this case the polarisation directions of photons rotate when propagating in the space 
[4] 1 . The rotation angle a, which characterises the CPT-violating signal, depends on the difference of between 
the source and the receiver of photons. It could be large as photons travel across the astronomical and cosmological 
distances, even though the coupling itself is very small. 

The cosmic CPT violation, if exists, could in principle be detected by the observation of the cosmic microwave 
background radiation (CMB) angular power spectra. Because the rotation of the polarisation vector of the CMB 
photons can covert part of the E-mode polarisation to B-mode polarisation, the non-zero TB and EB cross correlation, 
which vanishes in traditional CMB theory, can be induced [6]. The rotation angle can be probed by observations of 
the TB and EB power spectra [7], and the first constraint on the rotation angle using CMB observations was obtained 
in [8]. 

Recent studies along this line have demonstrated that current CMB experiments have the sensitivity to probe the 
rotation angle at the level of 0(l°)[9-22]. 

As quite a few B-mode CMB experiments are ongoing or being planned, we expect to have high quality observations 
of the TB and EB spectra, in addition to the TT, EE and TE spectra in the near future. Therefore we need a fast and 
accurate method to measure the rotation angle using these spectra. Traditional methods, e.g., [8-11, 13, 14, 17, 19, 20], 
which are based on the global fitting using the Markov Chain Monte Carlo (MCMC) method, are time consuming and 
computationally expensive. Alternative methods, which are more efficient, exist, e.g., [12, 16, 18], but are restricted 
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1 There are other mechanisms that can break the CPT symmetry, e.g., the circularly-polarised gravitational wave background (GWB) [5]. 
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to the special case where intrinsic BB vanishes. As the hint of a non-zero BB being discovered [24] and investigated 
[25], an efficient, and more general method including the intrinsic BB spectra are naturally needed. In this paper, we 
develop such a new method to probe the rotational angle based on a linear regression algorithm. We first validate the 
method using the simulated CMB data, before applying to recent CMB observations. 

The paper is outlined as follows. In section II, we shall present the methodology, including the tests using mock 
data. Section III shows the result of the measurement, followed by a discussion and conclusion section. 


II. METHODOLOGY 

In this section, we shall develop the methodology by first presenting the formulism, and then applying it to the 
simulated CMB data for a validation test. 


A. Formalism 


The rotation angle a(n) is anisotropic in general and can be decomposed into an isotropic part and an anisotropic 
one [26], namely, 


i(n) = a + 6 a (A) . 


(1) 


The CPT symmetry is violated by the background a, and the fluctuation <5a(n) will bring distortions to the CMB 
spectra, similar to the lensing effect. The technique based on the quadratic estimator to detect the anisotropies of 
a(n) has been developed in Refs. [27-29]. If the rotation angle itself has a spectrum, the rotated spectra of CMB can 
be calculated analytically [26]. 

The general relation between the unrotated and rotated angular power spectra Cp is [9, 26, 30], 
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and is the angular power spectrum of the rotation angle, say, 
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( 5 ) 


The overall constant A is related to C Q (0) via A = exp[—2C Q (0)]. Note that in the limit that the anisotropic part 
vanishes, i.e., da = 0, we have A = 1, f{t) = Cf E — Cf B and we get the familiar relation for a constant rotation 
angle [8, 10, 11, 13, 14, 17, 19, 20], 
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Prom Eq (2), it is straightforward to eliminate the unrotated C/s and obtain the following two linear relations 
between the rotated s, 


y^#TB,rot i /o — \ /'-yTE.rot /-» 

C £ — tan(2o;) C t ’ =0 

2C £ ’ — tan(4a) ( C i ’ — C £ ’ J = 0 (7) 

Note that these relations generally hold, even if the anisotropy part of the rotation angle exists (since it drops out), 
and they are cosmology-independent simply because the theoretical, unrotated Ct s do not show up here. This, in 
principle, allows the estimation of a directly from data. In what follows, we shall develop a new method based on the 
linear regression to estimate a, validate it using simulated data and apply it to the real CMB observations. 

The observed angular power spectra, C ° bs , are the rotated ones with measurement uncertainties. Therefore, the left 
hand side of Eq (7) do not strictly vanish when C\ ot is replaced by C£ bs , but with residues At and E^, namely, 

/\i = ’ — tan(2ai) C £ ’ 

= 2 - tan(4a) ( C £ - C £ J (8) 

It is a typical linear regression problem to find a from the above relation, and similar problems have been studied 
extensively in other fields of astronomy [31, 32]. Here we briefly review these methods. 

For a linear regression problem 2 , namely, 


Vi = @Xi, (9) 

where the ( Xi , yi) pair denotes the ith measurement with errors and the slope /3 is the free parameter to be determined. 
The developed linear regression methods include OLS (Ordinary Least Squares), BCES [33], FITEXY [34, 35], MLE 
(Maximum Likelihood Estimates) [36] and the sophisticated full Bayesian Gaussian mixing method [32]. Ref. [31] 
extensively compared the behaviour of all these methods based on the simulated data obeying various distributions, 
and concluded that the BCES, FITEXY, and the Bayesian methods work equally well. Here we adopt the FITEXY model 
due to its robustness and simplicity. 

The FITEXY method was originally proposed in [34] and modified by [35]. This method determines /3 by minimising 
the following y 2 (we ignore the correlations among data points for the moment), 
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Here (Tint accounts for the intrinsic scatter in the measurement of y , and its value is to be determined by iteration 
such that the reduced y 2 for the best fit model equals one, i.e., y 2 ed = Xbf/(-^ — 1) = 1- This is to guarantee that 
the dataset is never under-fitted in order to avoid the bias. As shown in Ref. [31], the OLS method (i.e., when the 
intrinsic scatter is not included) can lead to severely biased result. 

In general cases where the measurements are correlated, i.e., 

cov(ay, Xj) ^ 0; cov(y l ,y j ) ^ 0; cov(xi, yj) ± 0, (11) 

we generalise Eq (10) to, 

X 2 = (AH)(^ 5“)(f). (12, 

where 

F " s c = ( cil cii ) < 13 > 


We set the intercept to zero here since it does not apply for our concerned problem. 
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and 


(C A a)«/ = cov(A f ,A f /), 
(C H s)«/ = cov(H*,3*/), 
(Cas)h- = cov(A*,H*/), 
(C 3 a)«/ = cov(S^,Ap). 


( 14 ) 


Note that Cas and Csa are not symmetric simply because cov(A^,Sp) ^ cov(Ap,S^) but C is symmetric. These 
covariance matrices can be directly calculated given the full data covariance matrices, namely 3 , 


cov (Af, A?) 
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cov [(7™, C™] — tan(2a) {cov \Cj E , C 1 ™] + cov [Q™, C™] } 
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cov [2(7™ — tan(4a) (C EE - Cf B ) , 2(7 EB — tan(4a) (C EE - C BB )] 

4 cov [ Cf B , (7 eb ] 

—2 tan(4a) {cov [C EB , C EE ] — cov [C EB , C BB ] + cov [(7 EE , C EB ] — cov [C BB , C EB ]} 

+ tan 2 (4a) {cov [(7 EE , C BE ] — cov [(7 EE , C BB ] — cov [(7 BB , C EE ] + cov [C BB , (7 BB ] } 

cov [Cj B — tan(2a) Cj E , 2Cf B — tan(4a) (C EE — C BB )] 

2 cov [<7/ b , (7 eb ] - tan (4a) {cov [Q tb , Cf e ] - cov [< Cj B , C'J? 13 ] } - 2 tan (2a)cov [Cj®, Cf B ] 
+ tan (2a) tan (4a) {cov \Cj E , C EE ] — cov [Cj E , C BB ] } 


cov (3f, A#) = cov (Af, 3p)(^-B-^') (15) 

So given the binned measurements of Cj B and (7™, or (7 EE , C BB and (7 EB , the rotation angle can be determined 
without degeneracy with any other cosmological parameters, which is one of the advantages of this method. 


B. Tests using mock data 


To test the validity of this method, we generate mock CMB data following [9] and assume a Planck sensitivity 
summarised in Table 1 therein. The noise power spectra is then estimated as, 


Nl = (A T 6»FWHM, c ) 2 exp 
N^ c = (A P 6»FWHM,c) 2 exp 


1(1 + l)#FWHM,c 
8 ln2 

£(£ + l)flpWHM,c 
8 ln2 


(16) 


where T and P denote the ‘Temperature’ and ‘Polarisation’ respectively, and $fwhm,c is the Full width at half 
maximum (FWHM) of the angular resolution for a given frequency channel c. The combined noise from all channels 
is then, 
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3 We drop the superscript ‘obs’ for brevity. 
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FIG. 1: The validation test result using mock data. Each row of panels shows the result for a given fiducial model of rotation 
angle (top: a = 0°; middle: a = 1°; bottom: a = —2°). In each row, the left (middle) panel shows the linear regression result 
using mock data of TB and TE (2EB and EE-BB), where the data points with error bars are mock CMB data; the red (black) 
line shows the input (reconstructed) model. The blue error band shows the 68% CL uncertainty of the reconstruction. The 
panels in the right column show the one-dimensional posterior distribution of a derived from the linear regression method for 
three fiducial models. The red and black curves show the result using TE+TB and all the CMB data (including EE, BB and 
EB) respectively. 


The observed dimensionless Ct s (without the l{i + 1)/27t factor) are, 

Q bs = Cf eo + Nt (18) 

We randomly displace the data points by the corresponding 1-tr error. For this mock test, we choose a Planck best- 
fit cosmology with the tensor-to-scalar ratio r = 0.1, and with several fiducial values of the rotation angle, say, 
a = 0°, 1°, —2°. The mock data points are shown in the left and middle panels in Fig 1. One can easily identify the 
slope in this plot visually, which is useful for a quick consistency check. For example, if the slope in the TB-TE data 
has the opposite sign in the EE-BB-EB plot, it might suggest an inconsistency in the dataset itself. We then apply 
the FITEXY method to measure the slope between TB and TE, and between 2EB and (EE-BB) (we further drop the 
symbol of C(; for brevity). The measurements (best-fit value and 68% CL uncertainty) are over-plotted with data 
points in Fig 1. Given the measurement of the slope, we obtain the constraint on d, which is shown in the right panels 
of Fig 1 and in the top part of Table I. As shown, our method works very well in all three cases, namely, the input 
models are accurately reconstructed. 


Fiducial model for mock data 

TE+TB 

EE+BB+EB 

ALL 

Qfld = 0° 

-0.015° ±0.054° 

0.016° ± 0.020° 

0.013° ±0.018° 

Qfid = 1° 

0.984° ± 0.054° 

1.015° ±0.020° 

1.011° ±0.019° 

Qfid = -2° 

-2.020° ± 0.058° 

-1.981° ±0.020° 

-1.985° ±0.019° 


TABLE I: The measurement using the mock data (best fit value with 68% CL uncertainty) on the isotropic rotation angle. 
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III. RESULT 

In this section, we apply our method to recent CMB measurements, including ACTPol [23], BICEP1 [15], 
BOOMERanG 2003 (B03) [21], QUaD [16, 18], and we show the result in Figs. 2, 3 and in Table II. From Fig 
2, we can see that the the constraints from ACTPol and QUaD are consistent with zero rotation, but the slopes for 
the datasets of BICEP1 and B03 are nonzero at more than 68% confidence level. Our measurement shows, 

B03 : a = 3.7° ±6.8° (TE + TB); a = -13.1° ± 4.6° (EE + BB + EB); a = -8.1° ± 4.5° (ALL) 
BICEP1 : d =-3.6° ±1.2° (TE + TB); a = -2.6° ± 0.8° (EE + BB + EB); a = -2.9° ± 0.7° (ALL) (19) 

To be more conservative, we follow [9, 11, 37] to marginalise over another rotation angle ry with the following Gaussian 
prior to account for the possible systematics, namely, 


?7bo3 = -0.9° ± 0.7°; t?bicepi = 0.0° ± 1.3° (20) 

With the systematics included, the constraint is diluted to, 

B03 : a = -7.2° ±4.6° (ALL + sys.); BICEP 1 : a = -2.9° ± 1.5° (ALL + sys.) (21) 

which is consistent with the published result in [9]. The combined dataset gives, 

B03 + BICEP1: a = -3.3° ± 1.4°(ALL + sys.) (22) 

This shows a 2.3 a signal of a non-zero rotation angle. 



FIG. 2: The a measurement result using recent CMB polarisation observations. In each row, the left (middle) panel shows 
the linear regression result using the data of TB and TE (2EB and EE-BB), where the data points with error bars are CMB 
data; the red (black) line shows the input (reconstructed) model. The blue error band shows the 68% CL uncertainty of the 
reconstruction. The horizontal red dashed lines shows TB = EB = 0. The panels in the right column show the one-dimensional 
posterior distribution of a. The red dashed, blue dash-dotted and black solid curves show the result using TE+TB, EE+BB+EB 
and combined data respectively. The vertical dashed lines show a = 0 to guide eyes. 
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FIG. 3: Same as Fig 2 but for the most recent BICEP2 and BICEP2/Keck data. 


This signal can be visually identified in the BICEP1 and B03 panels of Fig 2. In both TB and EB panels of 
BICEP1, we can see a consistent negative slope at more than lcr level. However, the TB panel of B03 does not show 
any significant slope, but its EB panel shows a negative slope at 3cr level. The origin of this slope needs to be further 
investigated by future CMB experiments. 

Recently the BICEP2 collaboration published the new measurement of the CMB polarisation data [24], and a 
joint analysis between BICEP2 and Keck observation (BK) has been performed [25]. These are the most precision 
measurement of the CMB polarisation. But unfortunately these data cannot be used to constrain the rotation angle 
because they were already self-calibrated [38], so that any nonzero rotation angle, if exists, has been removed from 
the maps. But as a consistency test, we apply our method on these data, and show the result in Fig 3 and in the lower 
part of Table II. As shown, the constraint from the BICEP2 data is very consistent with a = 0, but the BK dataset 
shows a slight preference to a non-zero rotation angle, although the significance is low (1.7cr). This signal shows up 
when using the EE, BB and EB spectra, which gives a 1.7cr signal without using TE and TB. This might suggest that 
the self-calibration operation performed for the BK data is not as complete as that for the BICEP2 data. 

We show the correlation matrix, which is the rescaled covariance matrix C defined in Eq (14), for the best fit 
rotation angle using QUaD and BICEP2 data in Fig 4 4 . In each panel, the upper right and lower left blocks show 
the correlation coefficients between £ bins for the same kind of data (Css and Caa respectively), while the rest two 
blocks are for the correlation between £ bins for different kinds of data (Csa and Cas respectively). As shown, the 
correlation is only non-negligible between neighbouring bins for the same kind of data, or between different kinds of 
data but for the same bins, which is expected. 

This constraint on the rotation angle using our new, efficient method is largely consistent with those using the 
MCMC global fitting methods, as recently performed in [9]. Fig 5 shows the consistency, namely, all the measurements 
lay on the diagonal red line within lcr error bars. The two measurements using B03 and BICEP1 respectively deviate 
from a = 0 noticeably, which is discovered by both works. 


IV. CONCLUSION AND DISCUSSIONS 

Constraining the rotation angle of the polarisation vector of the CMB photons is an effective way to test the 
CPT symmetry on cosmological scales. Given that the CMB polarisation measurements are getting more and more 
accurate, and that several B-mode experiments are ongoing or being planned, it is timely to develop accurate, and 
efficient methods to probe the rotation angle. 


Note that the correlation matrix depends on a, so we need to specify the value of rotation angle for the illustration. 
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QUaD 


BICEP2 




FIG. 4: The correlation matrix for the best fit model using the QUaD and BICEP2 polarisation data respectively. The upper right 
and lower left blocks show the correlation coefficients between i bins for the same kind of data (Css and Caa respectively), 
while the rest two blocks are for the correlation between t bins for different kinds of data (Cea and Cas respectively). As 
shown, the correlation is only visible between neighbouring bins for the same kind of data, or between different kinds of data 
but for the same bins. 



FIG. 5: A visual comparison of the constraint on the rotation angle between this work and [9]. The red line illustrates the 
perfect consistency, and the horizontal and vertical dashed lines show a = 0. 


In this paper, we have developed a new method for this purpose based on the linear regression algorithm. Compared 
to previously established methods, our new method is computationally efficient (it takes several minutes on a laptop), 
accurate, intuitively transparent, and more generally applicable. This new method has successfully passed the vali¬ 
dation tests using mock CMB data before applied to CMB observations, including ACTPol, B03, BICEP1, QUaD, 
as well as BICEP2 and BK. We find that the B03 and BICEP1 samples prefer a nonzero rotation angle even when 
the systematics are included, and a combination of these two surveys gives a = —3.3° ± 1.4°, which is a 2.3cr signal. 
However, ACTPol and QUaD support a zero rotation angle. BICEP2 and BK provide the most stringent constraint 
on the rotation angle, namely, the uncertainty reaches ~ 0.1° level, but unfortunately they cannot be used to probe 
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Experiments 

TE+TB 

EE+BB+EB 

ALL (ALL+systematics) 

Published Result 

ACTPol [23] 
B03 [21] 
BICEP1[15] 
QUaD [16, 18] 

0.17° ±0.84° 

3.7° ±6.8° 

-3.6° ± 1.2° 

0.86° ± 1.79° 

0.34° ± 0.30° 

-13.1° ±4.6° 

-2.6° ±0.8° 

0.41° ±0.52° 

0.32° ±0.28° 

-8.1° ±4.5°(-7.2° ±4.6°) 
—2.9° ± 0.7°(—2.9° ± 1.5°) 
0.40° ±0.52° (0.40° ±0.72°) 

0.22° ± 0.32° [39] 
-4.63° ±4.16° [9] 
-2.69° ± 1.52° [9] 
0.59° ± 0.64° [9] 

BICEP2[24] 

BK[25] 

0.28° ± 0.27° 

0.09° ±0.19° 

0.08° ±0.17° 

0.17° ±0.10° 

0.11° ±0.17° 

0.15° ±0.09° 

0.12° ±0.16° [9] 


TABLE II: The measurement (best fit value with 68% CL uncertainty) on the isotropic rotation angle using different experi¬ 
ments. The numbers quoted in the parenthesis in the right column are those corrected with systematics. 


for the rotation angle since the datasets have been self-calibrated. 

With the high-quality polarisation data available in the near future, we would be able to constrain the anisotropic 
rotation angle 5a, as attempted by [9, 26-30]. Since 8a is degenerate with a using traditional methods, our method 
is efficient to break the degeneracy since a is not degenerate with any other cosmological parameters, including 5a, 
in our prescription. An extension of this work is to develop a new algorithm to probe the anisotropic rotation angle 
using the regression method and a principle component analysis, which is left for a future study. 
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